CLASSICAL COULOMB SYSTEMS : 
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Abstract 

From the laws of macroscopic electrostatics of conductors (in particular the existence 
of screening) taken for granted, one can deduce universal properties for the thermal fluctu- 
ations in a classical Coulomb system at equilibrium. The universality is especially apparent 
in the long-range correlations of the electrical potentials and fields. The charge fluctua- 
tions are derived from the field fluctuations. This is a convenient way for studying the 
surface charge fluctuations on a conductor with boundaries. Explicit results are given for 
simple geometries. The potentials and the fields have Gaussian fluctuations, except for a 
short-distance cutoff. 
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1. INTRODUCTION 



A salient property of matter is the screening effect : matter "in thermal equilib- 
rium does not tolerate any charge inhomogeneity over more than a few intermolecular 
distances". W In the present paper, we consider those systems which can be described 
as made of charged particles, interacting through Coulomb's law, to which equilibrium 
classical (i.e. non-quantum) statistical mechanics is applicable (for instance, electrolytes); 
then, screening has especially rich and simple consequences. From the existence of screen- 
ing, taken as granted, it is possible to deduce quantitative properties of the correlation 
functions. For instance, in the bulk, the charge-charge correlation function obeys the 
well-known Stillinger-Lovett sum rules; ^ other quantitative properties hold near walls or 
interfaces. ( 3 ) There is a delicate interplay between the statistical mechanics of correlations 
and the validity of macroscopic electrostatics. 

In the standard approaches, the focus is on the charges. One assumes that an external 
charge is screened and, using linear response theory, one obtains information about the 
charge correlations. The purpose of the present paper is to revisit the subject with the 
focus put on the electrical potential, the electrical field, and their correlations. From this 
alternative point of view, it is possible to rederive known results in a way that we believe 
to be often simpler and also to obtain new results. The present method , which relies 
on the validity of macroscopic electrostatics, is especially appropriate for studying surface 
properties such as surface charge correlations. 

The Coulomb systems considered in the present paper form a large class of models. 
The one-component plasma is a model such that identical particles of one sign (say positive) 
are immersed in a neutralizing inert continuous background of the opposite sign. The 
two-component plasma is made of two species of oppositely charged particles; a classical 
theory is possible only if some short-range interaction is also present, to prevent oppositely 
charged particles to collapse on each other. One could also consider more complicated 
models, with more than two species of particles, a non-uniform one-particle potential 
added, etc.. . We shall only require that the model be a conductor, in that sense that 
the laws of macroscopic electrostatics are assumed to be obeyed for length scales large 
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compared to the microscopic characteristic lengths of the model; for instance, an additional 
localized charge Q brought into the fluid should be perfectly screened, i.e. surround itself 
with a microscopic polarization cloud of charge —Q. 

We shall also consider models with a two-dimensional Coulomb interaction : in two 
dimensions, the potential at a distance r from a unit point charge is — ln(r/L) (where L 
is some fixed length) instead of the familiar 1/r which holds in three dimensions. Two- 
dimensional Coulomb systems made of oppositely charged particles exhibit the Kosterlitz- 
Thouless phase transition S 4 ^ 1 while they are conductors above some transition temperature, 
they become dielectrics below that temperature. A similar transition occurs for one- 
dimensional systems with a logarithmic interaction.^ The considerations of the present 
paper apply only to the conducting phase, i.e. above the transition temperature. 

The present approach aims to provide an easy to visualize physically reasonable picture 
of classical Coulomb fluids. No attempt is made towards mathematical rigor. Since a basic 
ingredient is macroscopic electrostatics, a number of words will be used with a macroscopic 
meaning. "Inside" or "outside" a Coulomb fluid will mean "at a distance from the walls 
large compared to the microscopic scale" . "Surface charge" will mean "charge in a layer 
of microscopic thickness under the surface". Etc... 

The two-point correlations inside and outside a three-dimensional Coulomb system are 
discussed in Sections 2 and 3, respectively. In Section 4, we establish a general method for 
deriving the surface charge correlations. Examples are given in Section 5. Section 6 is about 
conducting surfaces. Section 7 is about systems with the two-dimensional logarithmic 
Coulomb interaction. In Section 8, it is shown that potential and field fluctuations are 
Gaussian (except for microscopic distances). 

2. INSIDE A COULOMB SYSTEM 

In this Section, we consider some region inside a three-dimensional classical Coulomb 
fluid at equilibrium. This fluid is not necessarily homogeneous. 
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2.1. Potential and Field Correlations 

Let <f>(v) be the microscopic electrical potential at point r. We wish to rederive the 
asymptotic formula for the potential-potential correlation function^ 6 ) 

/?(0(r)0(r')) T ~ - 1 Ir-r'l^oo (2.1) 

where (3 is the inverse temperature and (. . ) T means a truncated equilibrium statistical 
average : {AB) T = (AB) - (A)(B). 

For deriving (2.1), we use screening and linear response. Let us put into the fluid an 
infinitesimal test charge q at r. The Hamiltonian Hq of the fluid must now be supplemented 
with the fluid-test charge interaction term which can be written as H' = q4>(r), where 
is the potential due to the fluid; it should be noted that <j> is defined as not including the 
potential due to q. By linear response theory, the average of this potential at some point 
r' is changed by 

{<P(r>)) q -( ( p(r')) = -Pq( ( p(r) ( p(r')) T (2.2) 

where < . . . > q means a statistical average computed with the full Hamiltonian Hq + H', 
while < . . . > is an average computed with H Q only. Our assumption about screening 
is that q surrounds itself with a polarisation cloud of microscopic size, of charge — q; the 
rest of the system is unchanged, except that, in the case of an insulated conductor, charge 
conservation requires that a charge +q spreads on the walls. Therefore, for |r' — r| large 
compared to the microscopic scale, the l.h.s. of (2.2) is the potential due to a point-like 
charge -g at r and a surface charge +q; at r' inside the conductor, that surface charge 
creates a constant potential q/C, where C is the capacitance. Altogether, for an insulated 
conductor, 

(«-(^')) = -^ + ^ (2-3) 
From (2.2) and (2.3), one obtains 

m^')) T = y 1 ^ - £ (2-4) 

For a grounded conductor, no surface charge appears, and the last term in (2.3) and (2.4) 
must be omitted. If the system becomes infinitely large, 1/C goes to zero. Thus, for 
infinite systems, one recovers (2.1) in all cases. 
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From our derivation, the "asymptotic" validity of (2.1) (for an infinite system) now 
has a more explicit meaning : (2.1) can be considered as an equality provided |r — r'| is 
large compared to the microscopic scale (the screening length). From now on, we shall 
write (2.1) and similar equations as equalities, with the understanding that there is some 
microscopic short-distance cutoff. This point of view, commonly used in field theory or 
in the theory of critical phenomena, is convenient for studying those properties which are 
independent of the microscopic detail, i.e. universal. 

Since the electrical field is E^(r) = — <9 M 0(r), the field-field correlation function is 
easily obtained from (2.4) (with or without the last term) as 

MF MF (t>W- ° 2 1 - 3(r-r%(r-r%-^(r-r') 2 (9 

while (E(r)) = since (0(r)) is a constant inside a conductor. 

2.2. Charge Correlations 

The charge density p(r) is related to the potential </>(r) by Poisson's equation A0 = 
—An p. Thus, by taking the Laplacian on r in both sides of (2.4), one obtains 

/3(p(r)0(r')) T = 5(r-r') (2.6) 

Of course, writing the r.h.s. of (2.6) as a delta function disregards some spreading and 
microscopic structure for small distances. But the simplified form (2.6) is enough for 
correctly giving, by integration upon r, 

J dr(p(r)<p(r')} T = 1 (2.7) 

which is the known Carnie-Chan sum rule.^ 7 -* 

In the case of an homogeneous fluid, the Carnie-Chan sum rule is equivalent to the 
Stillinger-Lovett^ 2 ) sum rule. We can retrieve that equivalence by taking the Laplacian on 
r' in both sides of (2.7), with the result 

-47r[3(p(r)p(r')) T = A5(r-r') (2.8) 

(Here the assumption that the fluid is homogeneous is necessary, because 5 actually stands 
for some peaked function of r - r'; if the fluid were no n- homogeneous, the microscopic 



width of that peak would be a function of r', and taking the Laplacian on r' would generate 
additional terms). From (2.8), one obtains 

j(p(r)p(r')) T dr' = (2.9a) 

and (after an integration by parts) 

~J <p(r)p(r')> T (r - r'fd(v - r') = 1 (2.96) 

Eqs (2.9a) and (2.9b) are the well-known Stillinger-Lovett sum rules; it is amusing that 
they can be written in the form (2.8), which, again, disregards the microscopic detail. 

3. OUTSIDE A COULOMB SYSTEM 

We now consider correlations involving at least one point outside the boundaries of a 
conductor. 

3.1. Potential and Field Correlations across the Boundary 

We want to compute the potential-potential correlation ((p(r)(j)(r')) T when r is inside 
the Coulomb fluid and r' outside. We can repeat the reasoning of Section 2.1, considering 
the response to an infinitesimal test charge q at r. The only difference is that the induced 
surface charge g, which appears in the case of an insulated conductor, creates at r' outside 
the conductor a potential qF(r') rather than q/C. This function F(r') is determined by 
macroscopic electrostatics : it is the potential created at r' when the conductor carries 
the total charge 1. Its explicit form can be obtained for simple shapes of the conductor. 2 
From (2.2), we now obtain 

/?<0(r)0(r')> T = - F(r') (3.1) 

2 In most textbooks about electrostatics, it is stated that the potential is continuous at 
the surface of a conductor. This would imply that (3.1) and (2.4) become identical when 
r' is on the surface. Actually, there is usually a potential difference across the surface of a 
conductor, due to the formation of an electrical double layer. However, a constant shift of 
the potential is of no effect on the electrical field. 
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The same result can be obtained by assuming that the test charge is put at r'; the total 
potential of the conductor then is qF(r') : because F can be regarded as an element of an 
inverse capacitance matrix, which is known to be symmetrical. 

The last term of (3.1) must be omitted in the case of a grounded conductor. It vanishes 
anyhow in the limit of an infinitely large conductor (such as, for instance, a conducting 
fluid filling a half-space) . 

From (3.1), one finds a field-field correlation which is again (2.3). 

3.2. Potential and Field Correlations Outside 

Finally, we want to compute the potential-potential corelation {<j)(r)(/)(r')) T when both 
points r and r' are outside the Coulomb fluid. We still consider the response to an infinites- 
imal test charge q at r. The total potential change at r' is some function gG(r, r'), with 
G(r, r') determined by macroscopic electrostatics : G(r, r ; ) is the total potential change 
at r' when a unit point charge is put at r. In the case of a finite conductor, for determin- 
ing G one must specify whether the conductor is insulated or grounded. G is explicitly 
computable for simple shapes of the conductor. 

The potential change due to the fluid only is q[G(r, r') — |r — r'| _1 ], and (2.2) gives 

/3(^(r)^(r')) T =^^-G'(r,rO (3.2) 
From (3.2), one finds for the field-field correlation 

P(E»(r)E„(r>)) T = ^ - G(r, r')] (3.3) 

4. SURFACE CHARGE CORRELATIONS 

In the present macroscopic approach, it is natural to introduce a surface charge den- 
sity cr, which will be associated to the electrical field discontinuity at the surface of the 
conductor. At some point r on the surface, 

4w(r)=£° ut (r)-4 n (r) (4.1.) 
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where the index n denotes the component normal to the surface (with the positive direc- 
tion defined as pointing towards the outside) and (r) is the limit of that field 
component as r is approached from the outside (inside). Therefore, the surface charge 
density correlation function is 

W(r)a(r')) T = J^([E™\r) - 4>)P™V) " E™(r>)]f (4.2) 

Using the field correlations (2.5) and (3.3), one finds the final result 

^ (r)CT(0)T = __^^|p| rir , esurface (4.3) 

Of course, this approach is valid only if the distance |r — r'| is large compared to 
the microscopic scale. Also, the surface charge density a has to be understood as being 
the microscopic volume charge density, integrated on some microscopic depth. Within 
these limitations, the computation of the correlation function (4.3) has been reduced to a 
problem in macroscopic electrostatics of conductors. 

A method for obtaining the charge correlations near the surface of a Coulomb system 
has been devised by Choquard et alS 8 ^ a few years ago. This method was presented as 
an approximation to the Debye-Hiickel approximation. However, it can be seen that the 
surface charge densities provided by the method of Choquard et al. have to be identi- 
cal to ours (their kernel G ex t is identical to our function G). Therefore, our approach 
rephrases the method of Choquard et a/., without any need for invoking the Debye-Hiickel 
approximation. 

5. A FEW SPECIFIC CASES 

In order to illustrate how the present method works, we now consider a few specific 
cases. 

5.1. Coulomb Fluid in a Half Space 

The surface of the Coulomb fluid is the plane xOy, which acts as an impenetrable wall, 
confining the particles to the half-space z < 0. The method of images gives, in the empty 
half-space, 

... .. 1 1 



where r* = (x,y,—z) is the image of r = (x,y,z). One readily finds from (4.3), on the 
plane xOy, 

^MO) T = - 87r2|r 1 _ r/|3 ( 5 -2) 

which is the result which had been obtained in ref. 9 in the equivalent microscopic lan- 
guage : the volume charge density correlation function (p(r)p(r')) T behaves asymptotically 
when (x — x') 2 + (y — y') 2 becomes large as F(z, z' , (x — x') 2 + (y — y') 2 ), where F is a 
function integrable upon z and z 1 obeying 

" L dz L dz ' F(z - z ' Ax - x ' f + (v - v ' f) = -*r\(v-m»* (B ' 3) 

5.2. Coulomb Fluid in a Ball 

A ball of macroscopic radius R, centered at the origin, is filled by a Coulomb fluid. 
The method of images now gives, outside the ball, 

1 R 1 R 

G(r, r' = I T \ ~ ~ H ~\ + — 5 - 4 

where the image of r has the coordinate r* = R 2 r/r 2 . The presence of the last term in (5.4) 
is appropriate for an insulated ball, i.e. if we want to study the fluctuations in an ensemble 
such that the total charge on the ball is constant, for instance the canonical ensemble; on 
the contrary, if the ball is assumed to be allowed to exchange charge with some reservoir, 
say at zero potential (grounded ball), the last term of (5.4) must be omitted. 
Using (5.4) in (4.3) gives, for an insulated ball, 

, Wr)ffM>T = _^ [ _L_ + J_ ] (5 . 5) 

where 9 is the angle between r and r'. In the case of a grounded ball, the last term in 
(5.5) must be suppressed, i.e. one must add the contribution +1/(4tt) 2 R 3 from the total 
charge fluctuations. 

This contribution from the fluctuations of the total charge Q can be alternatively 
derived by assuming that, for typical configurations, Q is uniformly spread on the surface. 
Then, the corresponding energy is Q 2 /2R. If one further assumes that the average of 
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this energy has the harmonic oscillator value (1/2)/? 1 , using a = Q/4kR 2 one finds the 
contribution 1/(4tc) 2 R 3 to (3 (a a). 

5.3. Coulomb Fluid in a Wedge 

Our method reproduces the results of Choquard et alS 8 \ These authors had expressed 
some doubts about the reliability of their result, because it was not obvious that their 
derivation, presented as an approximation to the Debye-Hiickel theory, was valid for a 
non-smooth surface such as a wedge. The present approach now puts things on the firmer 
basis of macroscopic electrostatics, that we believe to be true, even in the case of a wedge. 

The results in ref. 8 are different from the ones which had been obtained by Jancovici 
et alS 10 ^ 1 It is now apparent that the Section Wedge in ref. 10 is erroneous 3 , as well as its 
quotation in eqs. (10a) and (10b) of ref. 3. 

6. CONDUCTING SURFACES 

This Section is about two-dimensional Coulomb fluids, with an interaction potential 
which is the usual 1/r one. The surface charge correlations can be studied by a slightly 
modified version of the method which was used in Section 4. The field correlation function 
in three-dimensional space is still given by (3.3), and the surface charge correlations can be 
still expressed in terms of the correlations between the field discontinuities on the surface. 

In the case of a conducting plane, one finds 

/%(rMrOr = - 4 ^_ r/|3 (6.1) 

(This result could also have been derived by considering a slab*- 3 ' 10 -' in the zero thickness 
limit). It should be noted that the plane result (6.1) is not identical to the half-space result 
(5.2). This shows that, in the case of a three-dimensional Coulomb fluid with boundaries, 
there must be some coupling between the surface charge and the volume charge densities. 

3 The mistake in ref. 10 can be traced to the assumption following eq.(A.5); this assump- 
tion is just not true. 

10 



In the case of an insulated conducting spherical surface of radius R centered at the 
origin, one finds 

where 6 is the angle between r and r'. In the case of a grounded surface, the last term of 
(6.2) must be suppressed. Therefore, l/(4n) 2 R 3 is again the contribution from the total 
charge fluctuations, just as for a grounded ball filled with a three-dimensional Coulomb 
fluid, in agreement with the assumption that in this latter case the fluctuations of the total 
charge are localized on the surface only. The rest of (a(r)a(r')) T however is different for 
a spherical surface and for a ball filled with the fluid, as seen in (6.2) and (5.5); again, in 
the case of a ball, there is some coupling between the surface and the volume. 

7. TWO-DIMENSIONAL COULOMB INTERACTION 

The above method and results can be easily transposed to two-dimensional or one- 
dimensional systems with a two-dimensional Coulomb interaction : the potential at r 
created by a unit point charge at the origin is — ln(r/L), where L is some given length. 
These systems are interesting for a variety of reasons. In some cases they are exactly 
solvable models. Some of them appear in the theory of random matrices. 

One now obtains for the potential-potential correlations inside an insulated two- 
dimensional conductor 

/3(0(r)0(r')) T = -ln^ T ^-l (7.1.) 

where C is the capacitance; for instance, for a disk of radius R, C~ x = — ln(R/L). Now, 
if the system becomes infinite, C _1 — > oo and (7.1) has no well-behaved thermodynamic 
limit. ( n ) However, the corresponding field-field correlation function is given by 

while (-E M (r)) = 0, and these expressions remain finite for infinite systems. The analog of 
the surface charge density is now a charge per unit length that we still call cr(r). Following 
the same other steps as in Sections 2, 3, 4, one finds again that obtaining the "surface" 
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charge correlations < aa > for two-dimensional or one-dimensional systems with loga- 
rithmic interactions reduces (for macroscopic distances) to a macroscopic electrostatics 
problem : finding G(r, r'), the total potential at r' outside a conductor when a unit point 
charge has been put at r (also outside). For the explicit calculation of G, a specifically 
two-dimensional tool is provided by the theory of functions of a complex variable and 
conformal transformations. 

For a Coulomb fluid in a half plane (the fluid is assumed to fill the y < domain of 
the xyO plane), one finds 

^)^')) T = - 2^- xr <7 ' 3) 

which is the result which had been given in ref. 9 in microscopic language. Similarly, on 
an infinite conducting line, one finds 

(a result which had been previously derived by a different method^ 12 -'). 

For a Coulomb fluid in a disk of radius R centered at the origin, one finds 

P(a(r)a(r')) = - — -L (7.5) 
2iT z (2Rsm \y 

where 9 is the angle between r and r'. This result is valid both for an insulated and 
a grounded disk. In the case of a three-dimensional ball, grounding gave to the surface 
charge correlation function (5.5) an additional contribution coming from the fluctuations 
of the total charge. It is remarkable that such an additional term does not occur in two 
dimensions. The reason is that the energy associated to a total charge fluctuation Q on a 
disk is infinite, since increasing the charge Q by bringing in an additional charge 5Q from 
far away costs an energy Q5Q dr/r, where the integral diverges. Thus, fluctuations of 
the total charge cannot occur. 

As an illustration of the method of conformal transformations, let us consider a 
Coulomb fluid in a wedge. We use polar coordinates (r, 6) . The fluid is assumed to 
fill the two-dimensional wedge 7 < 9 < 2-7r; the domain < 9 < 7 is empty. For finding G 

12 



in the wedge < 9 < 7, we can start from its half-plane (7 = re) expression in terms of 
complex coordinates [z = rexp(z6*)] 

G 7=7r = ln||^| (7.6) 

z — z 

and make the conformal transformation z — > Z^ 7 which relates the half plane and the 
wedge. Thus, 

_ 1 ln r 2a + r' 2a - 2r Q r /Q cos a(0 + 9') 
where a = rr/'j, and one finds 

2 / /\a-l 

j3((T(r)ff(r')> = -t^t^— 1 rrr (7.8) 

where s = +1(— 1) if the two points r and r' are on the same side (different sides) of the 
wedge. 

Other specific cases are solved in refs. 13 and 14. 

8. GAUSSIAN BEHAVIOR OF THE POTENTIAL AND FIELD FLUCTU- 
ATIONS 

In Sections 2 and 3, we have obtained the two-point correlations for the electrical 
potential and field, except for their microscopic distance behavior. We now want to prove 
the stronger statement that the fluctuations of these quantities are "almost" Gaussian, with 
the above mentioned covariances (the meaning of "almost" will be explained hereafter). 

Showing that the fluctuations are Gaussian is a step towards better understanding 
why the finite-size corrections to the free energy of a Coulomb system^ 15 -* are similar to 
those of a critical field theory : the Gaussian model. 

The basic ingredient of the proof (already used in wave- number space in ref.15) is 
the assumption that the screening response (2.3) remains of the same form, linear in g, 
even if q is no longer infinitesimal but is an external point charge of arbitrary magnitude. 
Now, a linear response to a finite (non infinitesimal) perturbation characterizes a Gaussian 
distribution. As a pedagogical example, consider a particle in a one-dimensional potential 
V(x). A simple calculation shows that if, in thermal equilibrium, the average displacement 
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of the particle under an additional force q is proportional to q, then V(x) is a harmonic 
oscillator potential, and the probability density exp[— (3V{xj\ is Gaussian. 

For a proof in the present case of a Coulomb fluid, we study the cumulants, considering 
for instance an insulated fluid with three-dimensional Coulomb interactions. The argument 
is as follows. Place at n — 1 points ri, r 2 , r n _i (chosen inside or outside the fluid) point 
charges q\, q%, q n -\. Let r n be some point, inside or outside; if inside, r n is assumed 
to be at a distance large compared to the microscopic scale from those points among 
(i"i, r2, r n _i) which are inside. Since macroscopic electrostatics says that the potential- 
charge relations are linear, the response of (j)(r n ) is given by an obvious generalization of 
what has been derived in Sections 2 and 3 : 

n-1 

(<f>M){ qi} = (0(r n )> + t 8 - 1 ) 

i=i 

where ( ){ qi } means an average computed with the full Hamiltonian including the inter- 
action term Yl7=i <7i0( r i)) while ( ) means the average when there are no charges qi. The 
coefficients a» are those ones obtained in Sections 2 and 3 : 

a,i = — - 1 + — if both rj and r n are inside 

— r n | 6 

di = —- 1 + F(r out ) if one member (r ou ^) of the pair (r i5 r n ) 

| r i r n | 

is outside and the other one inside 
a,i = —-. 1 + G(ri, r n ) if both r$ and r n are outside (8.2) 

| r i r n | 

By definition, 

(0(rn)) te} = ^ J™ (8-3) 

Rewriting the r.h.s. of (8.3) as a logarithmic derivative and combining (8.1) and (8.3), one 
obtains 
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Pdq r . 

Since the cumulants (or truncated correlation functions) (. . .) T of the potential <j> are 
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defined by the expansion in powers of q\ 

m 1 \m 2 \ ... m n \ 

(8.5) 

(where the prime indicates the absence of the term with all m's simultaneously vanishing), 
a comparison between (8.4) and (8.5) shows that all cumulants of total order mi + m 2 + 
. . . + m n higher than 2 and with m n = 1 vanish. 

The above result can be rephrased as follows. Consider the cumulant (</>(ri)</>(r2) 
. . .4>{y p )) t ', with p > 2. That cumulant vanishes provided that one point at least, say r p , 
is at a distance large compared to the microscopic scale from all those other points which 
are inside the fluid; these other points however might be at arbitrary distances (including 
zero) from each other. Similar properties hold for the electrical fields. 

Were it not for the restriction about the position of r p , the vanishing of the cumulants 
of order larger than 2 would imply that the 4>(ri) are jointly Gaussian. With the restriction 
added, the statistics become "generalized Gaussian" . This restriction actually is not very 
drastic, or at least is of a kind which currently occurs in field theory. In the simple case of 
a homogeneous infinite system, the covariance /3(</>(ri)</>(r2)) = l/l r i — r 2| corresponds to a 
Gaussian probability density proportional to exp[— ((3/8tt) f (V</>) 2 <i 3 r], a form which is not 
unexpected since (l/87r)(V</>) 2 is just the Coulomb energy density expressed in terms of 
the electrical field — V0. It is well-known that this Gaussian field theory actually becomes 
singular for short distances, and the average ((/>(ri)</>(r 2 ) . . .(f)(r p )) is not properly defined 
when there are coincident points; some regularization is necessary. A most important 
property is that, at least in the absence of coincident points (i.e. when all distances 
are large compared to the microscopic cutoff), Wick's theorem can be straightforwardly 
used for expressing (0(ri)0(r 2 ) • • .(p(r p )) in terms of the 2-point function (0(ri)0(r 2 )) = 
i//?|n - r 2 j. 

The Gaussian behavior of the electrical fields follows. 

Our result about the higher-order cumulants had been previously obtained, by a dif- 
ferent method, in the special case of an infinite homogeneous one-component plasma/ 16 ) It 
should be noted that the quantities V^ ut which were shown to have Gaussian fluctuations 
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in ref. 6 are not the full potentials <fi considered here. 

The present considerations can be easily transcribed to the case of a two-dimensional 
(logarithmic) Coulomb potential. 

9. CONCLUSION 

The universal validity of macroscopic electrostatics implies in particular universal 
response properties for conductors. Whenever classical statistical mechanics is applicable, 
these response properties imply in turn universal properties for the equilibrium thermal 
fluctuations. Independently of the detail of the microscopic constitution of the conductor, 
for length scales large compared to the microscopic ones, the potential and field fluctuations 
are Gaussian with universal covariances, (outside the conductor, things depend on its 
geometry, but again not on its microscopic constitution). The surface charge fluctuations 
are also universal, depending only on the geometry. The volume charge fluctuation also 
have universal properties, but in a less explicit way, since the universality then appears 
only in the form of sum rules. 

Unfortunately, when quantum mechanics is used, the universal response properties 
no longer imply simple universal static (equal time) correlations, because of the more 
complicated form of the linear response theory. One can write quantum sum rules for the 
correlations in the special case of the one-component plasma* 17 ' 10 ), but these sum rules 
involve the non-universal plasma frequency. For more general quantum models, no sum 
rules are known for static or (real) time-dependent correlations. 
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